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1 Introduction 


A control system is said to be robust if it retains an assigned degree of stability or performance under 
perturbations. Robustness is considered as an elementary aspect in the analysis and design of control 
systems. There are several research lines for different models of uncertainties. Robust control under 
structured perturbations is an active area in recent yearst^Y[ll]_ 


The parametric robust stability problem was initiated by Kharitonov(1978)[^f who gave a sufficient 
and necessary condition for the robust stability of interval polynomials. Chapellat et al.(I989)[®l proved 
the Box Theorem for the multilinear model P{s) = U{s)V{s) + X(s)y(s). 


When considering multivariable systems, the robust stability problem becomes more complicated^" 5. 6. 7. 8. 9]_ 

Ref. [7] and Ref. [8] discussed the scalar matrix case and pointed out that a finite test does not hold. Ref. 

[5], [9] presented an n-dimensional result for polynomial matrices. Ref. [3] gave a 2n-dimensional result 
for the model B{s)A{s) + D{s)C{s), where B{s),D{s) are interval polynomial matrices and A(s),C'(s) 
are fixed ones. 


The purpose of this paper is to improve the robust stability criteria for multivariable interval control 
systems obtained by Kamal and Dehleh. By exploiting the uncertainty structures, we are able to reduce 
the computational burden in checking robust stability. 


2 Preliminaries 

Definition 2.1 Given an interval polynomial set r(s) = U £ its Kharitonov 

vertex set is ry(s) = {ri(s), r 2 (s), r 3 (s), r 4 (s)}, and its Kharitonov edge set is rE{s) = {Ari(s) + (1 — 
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>^)rj{s), (i,j) e {(1,2), (2,4), (4,3), (3,1)} A G [0,1]}, where 

ri (s) = Tq + rfs + + rf+ ... 

?’ 2 (s) = Tq + rf s + + rf + ... 

rsls) = Tq + rfs + r^s"^ + r]^+ ... 
r4(s) = + rYs + + rf + ... 

Definition 2.2 A polynomial is stable if all its zeroes belong to the open left half plane. A polynomial 
matrix is stable if its determinant is stable. 

Lemma 1 (Box Theoreml^l) Consider the polynomial family A(s) = {i5(s,p) = Fi(s)Pi(s) + 

... + Fm{s)Pm{s), Piis) are interval polynomials, Ai(s) are fixed polynomials, i = 1,... ,m} and sup¬ 
pose A(s) is without degree dropping. Then A(s) is Hurwitz stable if and only if A£;(s) is Hurwitz 
stable, where A£;(s) = ^i('S)ArPi(s) + Fiis)Ep^(s) + Fiis)Kp,is)} and Kp.{s), Ep^(s) 

are the Kharitonov vertex set, Kharitonov edge set of Pi{s), respectively. (Here we take X]i=r /* = 0 if 
r > t). 

Let A{s) = (ay (s))„xn = ('^b(s))„xn matrices, where A{s) = (aij(s))„x„ means 

that the matrix is an n x n one and its tj-th entry is aij{s). Given two interval polynomial matrices 
B{s) = {bijis))^^^,D{s) = (dy(s))„xn’ i'®'> bij (s), dij (s) are interval polynomials. Let 


B{ij) = [b,j{s) : b,j(s) = b[f + b[f s + .. 


(1) 

D{ij) = |dy(s) : d^j{s) = d[f + d[f s + .. 

. + dJV, d«G[4\dJ^]} 

(2) 

Bis) = {i?(s) : B{s) = 

(^*j('5))„xn ’ bijis) G i3(ij)| 

(3) 

Vis) = {Ois) : Dis) = 

((iy(s))^xn; SD(ij)| 

(4) 


Denote Kharitonov vertex sets and Kharitonov edge sets of D{ij) as Bviij),Pviij), PE{ij) and 

Vpitj) respectively. 


3 Main Results 


In this part, we consider the uncertain family 

A4(s) = Bis)A{s) + V{s)C{s) (5) 

Let Sn denote the set of all bijections of the set (1,..., n} onto itself. The following theorem was given 
by Kamal et al.(1996). 

Proposition 1® A4(s) is stable for all B(s) G B(s), D{s) G D(s) if and only if A4(s) is stable for all 
B{s) G BEis),D{s) G Ve{s), where 

( b,j{s) €BE{ij), j = (T{i) 'I 

Be{s) = < B{s) & B{s) : bij{s) e Bv{ij), j <j{i) > (6) 

[ a G J 

( d,j{s) €VE{ij), j = (j{i) 'I 

Pe{s) = < D{s) G V{s) : dy (s) G T>v{ij), j 7 ^ cr{i) ^ (7) 

[ a G J 

As a generalization of the Generalized Kharitonov Theorem, Proposition 1 addresses the robust sta¬ 
bility of A4(s) and reduce it to 4^" (n!)^ robust stability problems involving determinants with 2n pa¬ 
rameters. In this paper, we will show that, by making use of the uncertainty structure information, 
the original robust stability problem can be reduced to X]"=o(^n)^4^" (n!) robust stability problems 
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involving n parameters, where In the sequel, the notation X xY stands for the Cartesian 

product, which is the set of all ordered pairs {x,y), where x € X,y G Y. The following lemma is needed 
before we present our main result. 

Lemma 2 A{s), B{s),C{s) and D{s) are defined as before. Let 


N^\s) = 


3fci, ^2 G {1,..., n}; / = 1 ,... , n 

fV«(s) = <1 B{s)A{s) + D{s)C{s) : 6,^, G S£;(^fcl); 6*; G Bv{il), I h 

diki G VE{ik 2 )]dii G Vv{il), I 7^ 

3fc G {1,..., n}; I = 1,... , n 
i3(s)2l(s) + D{s)C{s) : biu x dik G {BE{ik) x Vv{ik)) U {Bv{ik) x VEiik)) 

bii X dii G {Bv{il) X 'Dv{il)),l ^ k 


( 8 ) 

(9) 


Then, for alH G {1,..., n}, N''^\s) is stable if and only if N^\s) is stable . 

Proof (Sufficiency) Suppose that N^\s) is stable , our aim is to prove that N^^\s) is stable. For any 
N{s) G fV(*^(s), there exist fci, ^2 G {1,... ,n} such that 


biki G Be (^^ 1 ) 5 bii G B\^ (^ 0 ? I ^ k\j I — l,..., 7 r 
diki G BEij'k ^), dii G By (^ 0 ^ ^ 7 ^ ^ 25 ^ — 1 ,..., u 


( 10 ) 


and 


iV(s) = 


/E. (6ifc(s)afci(s) + difc(s)cfci(s)) 

^(0;^ (^zfc ('^)afci (s) T difc (s)c;^l (s)) 

V Efc (^nfc ('^)n/cl (s) + dji/,;(s)Cfci (s)) 


Efe (^lfc('S)afen(s) +'^lfc(s)Cfcn(s)) \ 
Efe ibik{s)aknis) + d^k{s)ckn{s)) 

Efc (^)nfcn (^) T djiki.^'j^kni^)) J 


( 11 ) 


By definition, N{s) is stable if and only if detA^(s) is stable. By Laplace formula, expanding the deter¬ 
minant along the t-th row, we have 


detA^(s) — ^ ) (^ifc ('^)afcl (s) -t- dik(^S^Ckl{s^'j Xu -t- . . . T ^ ) (biki^^^^kni^^ T diki^S^Ckni^)^ 
k k 


where Nij is the algebraic complementary minor of the ij-th entry of N{s), which is independent of 
bii{s ),..., bin{s), dii{s), ..., din{s). By a simple manipulation, we have 

detA^(s) = b^Jis) {YTk=i aok{s)N,k) + EJ=i d^J{s) {YTk=i Cok{s)N,k) 

= biki{s)6i{s) + dik^{s)52{s) + ( 53 ( s ) 


where (5i(s) is a term which is independent of 6ii(s),..., fcm(s), dii{s ),..., din(s). Denote 

( bik^{s) X dik^is) G {BE{iki) x Bv{ik 2 )) U {Bv{iki) x ^^(ife)) 

< bii{s) G Bviil),l^ki,l = l,...,n (12) 

[ dii{s) &Bv{il),l k 2 ,l = 1, ■ ■ ■ ,n 

By Lemma 1, detA^(s) is stable for N{s) satisfying (10),(11) if and only if detA^(s) is stable for N{s) 
satisfying (11),(12). Since N^\s) is stable , it shows that N{s) is stable. Hence, sufficiency is proved. 
(Necessary) This part is obvious, since N^\s) C 


Before presenting our main result, we first define some notations. Let T„ be the set of all functions 
from the set {1,..., n} into itself. Take 

( bij{s) X dij{s) G {BE{ij) X Bv{ij)) U {Bv{ij) x BE{ij)) j = vii) 

BB = < B{s) X D{s) : bij{s) x dij{s) e Bv{ij) x Bv{ij) j 

I ?? G T„ 
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and 


BeVe = {B{s) X D{s) : B{s) x D{s) S Be{s) x Ve{s)} . 

By Proposition 1 and Lemma 2, we get the main result of this paper. 

Theorem 1 M{s) is stable for all B{s) G B{s),D{s) G I?(s) if and only if A^(s) is stable for all 
B(s) X D(s) G BT) n BeBe- 

Proof By Proposition 1, M{s) is stable for all B{s) G B{s),D{s) G T>{s) if and only if M{s) is stable 
for all B{s) x D{s) G BeBe- For any n(s) G M{s) with B{s) x D{s) G BeBe, fixing the rows from the 
second to the n-th and applying Lemma 2 to the first row, we have 

A^(s) is stable for all B{s) G B{s), D{s) G T>{s) if and only if 

M{s) is stable for all B{s) x D{s) G BT>i n BeBe 

where BBi C B{s) x B{s) satisfies (12) when i = 1. Repeating the same procedure for the remaining 
rows consecutively, we get the conclusions. 

Remark 1 Theorem 1 shows that, the criterion for the robust stability of interval polynomial matrices 
can be greatly simplified. The dimension of the minimal testing set is reduced to n. 

Remark 2 Theorem 1 establishes the minimal dimensional testing criteria for robust stability of 
multivariable control systems. By making use of the uncertainty information of B{s) and B{s), i.e., not 
only making use of the uncertainty information in B{s) and B{s) separately (this is exactly what Kamal 
and Dehleh did), but also making use of the uncertainty information collectively, we are able to improve 
the robust stability results given by Kamal and Dehleh. 

Remark 3 Theorem 1 allows us to improve the main results of Kamal and Dehleh. 


4 Applications 


Based on Proposition 1, Kamal et al.(1996) address several robustness problems as follows. 

Proposition 2([^1) Let P{S) = B{S)D~^{s) and K{s) = yl(s)C'“^(s) be right coprime factorization 
and left coprime factorization of P{s) and K{s) respectively over the ring of polynomial matrices. Let 
B{s) G B{s) and D{s) G B{s), and let K{s) be fixed. Then the closed loop system in Fig. 1 is stable for 
all B{s) G B{s),D{s) G B{s) if and only if it is stable for all B{s) G BeIs), D{s) G Be{s). 



Fig 1: Feedback Systems 


Let H^'^{C+) be the space of matrix valued functions F{s) that are analytic in €+ and bounded on 
the jw-axis with the norm IjF’Hoo = sup^^g,^ (Tniax(F’(jw)). 

Proposition 3([^1) Let Q(s) C iL^^”(C+) be a family of proper interval transfer function matrices. 
Let G(s) = B(s)D~^(s) be a right coprime description of G(s) over the ring of polynomial matrices, with 
B(s) column reduced; then 

1) l|G||oo < 1 for all B(s) G B(s), D{s) G B{s) if and only if ||G||oo < 1 for all B{s) G Be{s), D{s) G Be{s). 

2) G(s) is strictly positive real for all B(s) G Bis),D(s) G B(s) if and only if G(s) is strictly positive 
real for all B{s) G BEis),D{s) G Pis(s). 
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A memoryless nonlinearity : [0, oo) x 7^" ^ 7?.” is said to satisfy a sector condition if for all y G 7?." 
and t > 0 , 

[As{t, y) - KiyflAsit, y) - K 2 y] > 0 (13) 

for some real matrices Ki and K 2 , where K = K 2 — Ki is a symmetric positive definition matrix. 

Proposition 4([^1) G(s) = B{s)D~^{s) is a family of interval proper stable systems, where B(s) G 
B{s), D{s) G 27(s), and is a set of memoryless nonlinearities that satisfy the sector condition globally 
with K > 0, then the system is globally stable if there is 77 > 0 such that G{oo) = 0 and 

Z^juj) =/+(!+ jujrj)KG{juj) > 0, Vw G TZ 
for all B{s) G Be{s),D{s) G VEis). 


By resort to Theorem 1, all testing sets above can be reduced to lower dimensional ones. 

Theorem 2 Under the same assumptions as in Proposition 2, the closed loop system in Fig. 1 is 
stable for all B{s) G B{s), D{s) G B{s) if and only if it is stable for all B{s) x D{s) G BV n BeBe- 

Theorem 3 Under the same assumptions as in Proposition 3, we have 

1) l|G||oo < 1 for all B{s) G B{s), D{s) G T>{s) if and only if ||G||oo < 1 for all B{s) x D{s) G BVCiBeBe- 

2) G(s) is strictly positive real for all B(s) G B(s),D(s) G T>(s) if and only if G(s) is strictly positive 
real for all B{s) xD{s)€BVn BeVe- 

Theorem 4 Under the same assumptions as in Proposition 4, the system is globally stable if there is 
?7 > 0 such that G{oo) = 0 and 

Z^juj) =/ + (!+ juj'q)KG{juj) > 0, Vw G TZ 

for all B{s) x D{s) eBVr\ BeVe- 


5 Example 


Recall the two-link planar manipulator considered by Kamal and Dehleh. The characteristic polynomial 
of the closed loop system is 

T(e) = det {M{9d)s^ + KdS^ + KpS + Kr) (14) 

where 9d G [0,7r/2] is the joint angle and 


Mi9d) 

Kd = 

Kp = 

Kr = I 


/ 3-1-2 cos{9d) 
1 -I- cos{9d) 


kdii 

kdi2 

kd2i 

kd22 

kpii 

kpi2 

kp2i 

kp22 

krll 

krl2 

kr21 

kr22 


1 -I- cos{9d) ^ 
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kdii e [6.07-6.07e, 6.07+6.07e] 
kdi 2 e [2.22 - 2.22e, 2.22 + 2.22e] 
kd 2 i G [2.22 - 2.22e, 2.22 + 2.22e] 
kd22 G [1.62 - 1.62e, 1.62 + 1.62e] 
fcpii G [6.12 - 6.12e, 6.12 + 6.12e] 
fcpi2 G [2.24 - 2.24e, 2.24 + 2.24e] 
kp2i G [2.24 - 2.24e, 2.24 + 2.24e] 
kp22 G [1.64 - 1.64e, 1.64 + 1.64e] 
krii G [5.11 - 5.11e,5.11 + 5.11e] 
kri 2 G [1.87 - 1.87e, 1.87 + 1.87e] 
kr 2 i G [1.87 - 1.87e, 1.87 + 1.87e] 
kr 22 G [1.37 - 1.37e, 1.37 + 1.37e]. 


Since for any fixed 9d G [0,7r/2], the inertia matrix 


M{9d) = 


1 


(-1 + cos(6»d)) 1 

In what follows, we take the uncertain inertia matrix as 


(1 + cos{9d)) ^ 1 0 

2 1 


M{s) = B{s)A{s) 

where A(s) = 2 1 )’ ^ ^ {ibij{s)) 2 x 2 ■ bij{s) G B{ij)} and 

B{ll) = s^ B(12) = 6i2s3, 6i2G[1,2] 

B{21) = 62 lS^ 521 G [-1,0] B{22) = 


( 15 ) 


Let V{s) = KdS^ + KpS^ + Kr, then 

V{s) = {{dij{s)) : dij(s) G D{ij)} , 

where D{ij) = {(ijj(s) = kdijS^ + kpijS + krij : kdij, kpij and krij are given in (15)}. Thus, the charac¬ 
teristic polynomial of the uncertain system is 


r(e) = det {B{s)A{s) + D{s )), B{s) G B{s), D{s) G V{s) 


(16) 


where 2 l(s) is a fixed matrix and B{s), D{s) are interval polynomial matrices. 


The Kharitonov polynomial vertex sets associated with B{ij) and D{ij) are as follows. 


= B{n), Bv{12) = {s3, 2 s3}, ^y(21) = {-s^ 0}, Bv{22) = {s^} = 5(22), 

+ ^rij^ + ^pij^ + ^rij^ \ 

k)f = I 


Vv{ij) = 


-I- s -I- 


where the upper bound and the lower bound of Kdij are denoted as . And the Kharitonov edge 

sets associated with B{ij) and D{ij) are as follows. 


5^(11) = {s3} = 5(11), 
5^(21) = {-As3; AG [0,1]}, 


BEiij) = 


Be{12) = {As3 + (1 - A) • 2s3 : a G [0,1]}, 
5b( 22) = {s3} = (22) 

Afc« + (1-A)fcg\ } 

Afc|) + (l_A)fc«, 

A^) + (1 - A)4\ 

Afc® + (1-A)4\ 
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By Proposition 1 given by Kamal and Dehleh, for any e, T{e) is robustly stable if and only if T(e) is 
stable for all B{s) x D{s) G Be{s) x VEis), where 

Be{s) xVe{s) = Ui\JU2 
( B{ll) Be{12) \ ( Dv(ll) De(12) \ 

I 5^(21) B{22) j \ De{21) Dv{ 22) j . 

( 5(11) Be{12) \ ( De{11) Dv(l2) \ 

V Be(21) 5(22) ) \ Dv(21) 5^(22) ) 

So, the robust stability of T(e) is reduced to 2 • 4® robust stability problems involving determinants with 
four parameters. In the case, the testing set is a four-dimensional set. In the sequel, we apply Theorem 
1 to simplify this problem further. That is to say, T(e) is robustly stable if and only if T(e) is stable for 
all 5(s) X D{s) G {BV) n BeBe, where 


bdUV^UVsUVljUVsUVgU Vy 


{BV) n BeVe 

( 5(11) 5£;(12) 

I 5£;(21) 5(22) 

/ ^( 11 ) Be{12) 

2 I 5y(21) 5(22) 

( 5(11) 5£;(12) 

^ I 5v(21) 5(22) 

( 5(11) 5v(12) 

I 5£;(21) 5(22) 

/ 5(11) 5v(12) 

I 5^(21) 5(22) 

/ 5(11) 5y(12) 

® I 5y(21) 5(22) 

/ 5(11) 5y(12) 

^ 5y(21) 5(22) 


\ / 5y(ll) 5y(12) 

J { Dv(21) Dv{22) 

\ ( De{11) Dv{l2) 

\ Dv{2l) Dv{22) 
\ ( Dv{ll) Dv{l2) 

j \ De{21) Dv{22) 

\ ( Dv{ll) De{12) 

\ Dv{2l) Dv{22) 
\ ( Dv{ll) Dv{l2) 

j I 5y(21) 5£;(22) 

\ ( De{11) Dv{l2) 

j I 5y(21) De{22) 

\ ( Dv{ll) De{12) 

) ^ \ De{21) Dv{22) 


Clearly, the robust stability of T(e) is reduced to 7 • 4® robust stability problems involving determinants 
with two parameters. In this case, the minimal testing set is two-dimensional. 


Remark 4: A test with a parameter could involve infinite tests without parameter, such as an edge 
means infinite points. The significance of our result is that it is the minimal dimension test for such 
problem as (5). Furthermore, we consider the uncertainty in inertia matrix M{0d) in Example, which 
did not take into consideration in Ref. [3]. 


6 Conclusions 


In this paper, interval transfer function matrices are studied. For a family of multivariable control 
systems having interval transfer functions matrices, a lower dimensional robust stability criterion has 
been established. This result improves the previous criterion obtained by Kamal and Dahleh, and can be 
used to establish improved lower dimensional robust performance criteria as well. 
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